Gluon saturation and baryon stopping in the SPS, RHIC, and LHC energy regions 
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A new geometrical scaling method with gluon saturation rapidity limit is proposed to study the 
gluon saturation feature of central rapidity region of relativistic nuclear collisions. The net-baryon 
number is essentially transported by valence quarks that probe the saturation regime in the target 
by multiple scattering. We take advantage of the gluon saturation model with geometric scaling 
of rapidity limit to investigate the net baryon distributions, nuclear stopping power and gluon 
saturation feature in the SPS and RHIC energy regions. Predications for the net-baryon rapidity 
distributions, the mean rapidity loss and gluon saturation feature in central Pb -I- Pb collisions at 
Life are made in this paper. 
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I. INTRODUCTION 

During the relativistic heavy-ion collisions, the fast va- 
lence quarks in one nucleus scatter in the other nucleus 
by exchanging soft gluons, leading to their redistribution 
in rapidity space. The net-baryon number is essentially 
transported by valence quarks that probe the saturation 
regime in the target by multiple scattering. 

Experimental heavy-ion investigations at the Large 
Hadron Collider(LHC) pay much more attention to the 
mid-rapidity region since ALICE covers rapidity up to 
1 2/1 = 2. It provides measurements of lower x than before 
at energy scale that is high enough to provide crucial 
tests of gluon saturation. Therefore LHC will provide 
more opportunities to study the gluon saturation feature 
at small x. 

At very high energies or small values of Bjorken vari- 
able X, the density of partons per unit transverse area, in 
nucleon or nucleus becomes so large that it would lead to 
a gluon saturation of partonic distribution. The existence 
of this phenomenon was confirmed in the experiments at 
HERA [l|, Q • The typical results from the experiments 
contain two parts: the small x problem and the geometric 
scaling. It was predicted by an effective theory, the CGC 
(Color Glass Condensate), which describes the behavior 
of the small x components of the hadronic wave function 
in QCD. In this paper, we use this CGC theory to discuss 
the questions of nuclear stopping and gluon saturation in 
relativistic heavy- ion collisions at SPS, RHIC and LHC. 

The kinetic energy of the relativistic heavy-ion colli- 
sions that is removed from the beam and which is avail- 
able for the production of a state such as the QGP (quark 
gluon plasma) depends on the amount of stopping be- 
tween the colliding ions. The stopping can be estimated 
from the rapidity loss experienced by the baryons in the 
colliding nuclei [3|-|9[ . If the incoming beam baryons have 
rapidity, yb relative to the CM(which has y ~ 0) and av- 
erage rapidity 



<y > 

the average rapidity loss is 



Vb rVb 

ydydN/dy/ / dydN/dy (1) 
Ja 



< Sy >= yb- <y> (2) 

Here dN/dy denotes the number of net-baryons (num- 
ber of baryons minus the number of anti-baryons) per 
unit of rapidity. The studies of net baryon distributions 
and nuclear stopping have been discussed by non-uniform 
collective flow model [I0i-tl3j . We will use gluon satura- 
tion to study net-baryon in this paper. 

Here we should mention a novel gluon saturation 
model proposed by Yacine Mehtar-Tani and Georg 
Wolschin [T3, EBl- An analytical scaling law is derived 
within the color glass condensate framework based on 
small-cou pling QCD in this model. Inspired by this 
model UAUM, we study the net baryon distributions of 
central collisions in the SPS and RHIC energy regions by 
introducing the effective quark mass and rapidity limit 
of gluon saturation region. Here we use a new geometric 
scaling method to define the rapidity limit of gluon satu- 
ration region. The important difference is the definition 
of geometric scaling between our model and Ref.[f4,15]. 
It realized [14, I5]that the geometric scaling is mainly 
about the peak positions of net-baryon rapidity distri- 
bution. But our model realizes the geometric scaling 
of net-baryon rapidity distribution is mainly about the 
gluon saturation rapidity limit. The net-baryon rapidity 
distributions and the mean rapidity loss in central Pb -I- 
Pb collisions at LHC are predicted in this paper. 



II. GLUON SATURATION MODEL WITH 
GEOMETRIC SCALING 

The ideas for the color glass condensate are motivated 
by HERA data on the gluon distribution function 
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The gluon density, xG{x,Q'^), rises rapidly as a func- 
tion of decreasing fractional momentum, a;, or increasing 
resolution Q. This rise in the gluon density ultimately 
owes its origin to the non-Abelian nature of QCD and 
that the gluons carry color charge. At higher and higher 
energies, smaller x and larger Q become kinematics ac- 
cessible. The rapid rise with Injl/x) was expected in a 
variety of theoretical works [14| - |23j . Due to the intrin- 
sic non-linearity of QCD, gluon showers generate more 
gluon showers producing an exponential avalanche to- 
ward small X. The physical consequence of this exponen- 
tial growth is that the density of gluons per unit area 
per unit rapidity of any hadrons including nuclei must 
increase rapidly as x decreases [H, [l^ . 

The net-baryon number of relativistic heavy ion col- 
lisions is essentially transported by valence quarks that 
probe the saturation regime in the target by multiple 
scatterings. During the relativistic heavy ion collisions, 
the fast valence quarks in one nucleus scatter in the other 
nucleus by exchanging soft gluons, leading to their re- 
distribution in rapidity space. The valence quark par- 
ton distribution at large x is well known, which corre- 
sponds to the forward and backward rapidity region, to 
access the gluon distribution at small x in the target nu- 
cleus. Therefore, this picture provides a clean probe of 
the gluon distribution, ip[x,pT) , at small x in the satu- 
ration regime. 

For symmetric heavy ion collisions, the contribution of 
the fragmentation of the valence quarks in the projection 
is given by the simple formula for the rapidity distribu- 
tion of interactions with gluon in the target EBl 



dN 
dy 



C 



2 

PT 



Xiqy{xi)ip{x2,PT) 



(3) 



here xiqu{xi) is the valence quark distribution of a nu- 
cleus, (p(x2,pt) is the gluon distribution of another nu- 
cleus, pt is the transverse momentum of the produced 
quark and y its rapidity. One important prediction of the 
gluon saturation with geometric scaling is the geometric 
scaling: the gluon distribution depends on x and pr only 
through the scaling variable Pt/QH^), here Q^{x) — 
A^/^Qpa;"^. Geometric scaling was confirmed experimen- 
tally at HERA [l|.The fit value A = 0.2 - 0.3 agrees 
with theoretical estimates based on next-to-leading order 
Balitskii-Fadin-Kuraev-Lipatov (BFKL) results [24 26]. 

The longitudinal momentum fractions carried, respec- 
tively, by the valence quark in the projectile and the soft 
gluon in the target are x = xi — [pr / ^/s)exp{y) , and 
X2 = (pt / \fs)exp(—y) . Then we can deduce the relation 
as follows VU, 



-2v 2 2 — 2j/ 

X = xit X2 = xe " , prp — X se " 



(4) 



The gluon distribution is defined as 



fix2,PT) = 



Pt 
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When we use the variables as Eq.4, we can take the 
relation as follows 

pI ^ {x-^s-e-yf ^ x^+^ _s_ J_ 

Q2(X2) v41/3Q2^-> e2(l+A)y ■ Ql ■ ^1/3 ^ ^ 



If we make ^^^7 ^ = x^^^ • e'^from Eq.6 we can take 
1 s 1 



(7) 



g2(l-HA)y q2 ^1/3 

and a geometrical scaling with rapidity is introduced as 



T ^ ln{^) - In A^/^ - 2{1 + X)y (8) 



Thus Eq.3 is given as follows 



dN 
dy 



C dx 
2^ 



xqv{x)ip{x'^^^ ,e'') 



(9) 



As mentioned before, the dependence of t on y in 
Ref. [14,15] is related to the peak position, but here we 
use the dependence of r on y to define the rapidity limit 
of gluon saturation. A rapidity variable is introduced in 
Eq.8 to discuss the gluon saturation rapidity region. The 
rapidity limit variable y is 

,2 



y = 



here yi, « ln{y/s/mn) is the beam rapidity with nucleon 
mass m„. It usually defines x ^ 0.01 as a small x region, 
and in this region, it is taken as the gluon saturation 
region. By taking r = ln{l/x)\x=o.oi : we can figure out 
the rapidity region of gluon saturation as follows 



y saturation 

= j^iyt-lnA'^') + 

The gluon distribution is 
I const 



2(1 + A) ^ Ql 



1 /J I ^ 

-{In-TTT - t\x=o.oi, 



(11) 



{x < 0.01) 



v{x2,Pt) 
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• (1 - X2)'^ {x > 0.01) 
(12) 



The valence quarks distribution is 



XQv (x) 



cx X 



,0.5 



{x < 0.01) 

oc (xUy + xdjj) (a; > 0.01) 



(13) 



So the net-baryon distribution which originates from the 
projectile is 

'cx exp{^^y} [0 < y < ysaturatton) 



dN 
dy 



X2)'^ {y > y saturation ) 
(14) 
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The contribution of valence quarks in the other beam 
nucleus is added incoherently by changing y ^ —y. The 
total rapidity distributions of the symmetry interaction 
systems are the summation of the contributions from the 
projectile and target, respectively 



dN , dN , , dN , , 

— U., = —{y) + —{-y) 



(15) 



Compared with Ref. [14,15], we will introduce effec- 
tive dynamic quark mass in this paper by substituting 
the transverse momentum pT as ^ (xy/se^y)"^ — m? , i.e. 

Pt — \/ {xy/se~yY ~ here m is the dynamic quark 
mass, which was given by baryon mass in Ref. [14,15]. 
The detailed study of dynamic mass with collision ener- 
gies will be given in the next section. 



III. THE NET BARYON DISTRIBUTIONS IN 
THE SPS TO LHC ENERGY REGIONS 

The integrated net-proton rapidity distributions are 
scaled by a factor of 2.05 |5| to obtain the net-baryon 
distributions. The rapidity distributions of net baryon 
for different energies of relativistic heavy nuclear colli- 
sions (Pb -I- Pb and Au + An) at SPS and RHIC are 
given in Figure 1. 

The estimated numbers of participants are 390, 315, 
and 357 for ,JJWn = 17.3, 62.4, and 200 GeV, respec- 
tively. The solid circles correspond to the experimental 
result of central collisions , and the real lines are the 
calculated results from the gluon saturation model with 
geometric scaling. It is found that our model describes 
the experimental data of net-baryon distributions very 
well when we discuss Pb-Pb center collisions at the SPS 
energy region and Au-Au center collisions at the RHIC 
energy region. The A = 0.2 and Qq = 0.05 GeV^ are 
fixed in our whole calculations. 

We show in Fig. 2 the computation resulting from our 
discussions, the centrality dependence of the rapidity dis- 
tribution at y^SNN = 200 GeV. The estimated numbers 
of participants are 280, 200, 114 and 54 for centralities of 
0% - 10%, 10% - 20%, 20% - 40%, and 40% - 60%(from 
top to bottom) in Au-t-Au collisions at RHIC energies of 
y/sNN = 200 GeV, respectively. 

By studying the experimental results with our model, 
we may get the conclusion for the net-baryon distribu- 
tion from SPS to LHC as follows: IThe gluon saturation 
model with geometric scaling may be a good theory to 
study net-baryon distribution at central collisions at SPS 
and RHIC. The contributions of spectator nucleons can 
be neglected when considering only the central collisions. 
The gluon saturation feature of central rapidity can be 
studied from our discussion. The values of central rapid- 
ity of gluon saturation (ysaturation) are 1.06, 2.013 and 
3.10 from Eq.ll for = 17.3, 62.4, and 200 GeV, 

respectively. 
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FIG. 1: Net-Baryon distributions of central collisions at SPS 
^/sNM = 17.3 GeV of Pb-Pb interactions and at RHIC 
ysjviv = 62.4 and 200 GeV of Au-Au interactions; the ex- 
perimental results are from [3-9]. 



(2) The net baryon rapidity distribution in central 
Pb-I-Pb collisions at LHC energies of ^/snn = 5.52 TeV 
is predicted by gluon saturation model with geometric 
scaling. The theoretical distribution is shown in Fig. 3 
for ysaturation = 5.86. The gluon saturation region is 
larger than those of RHIC and SPS. It is found that the 
separation of two symmetric peaks of net-baryon is much 
wider than that of SPS and RHIC. 

(3) The mean rapidity loss < Sy >— yb— < y > is 
shown in Fig. 4. We show that the dependence of mean 
rapidity loss increase on as 



< Sy >= 1.548 •ln(7;h) -f 0.036 



(16) 
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FIG. 2: Rapidity distributions of net baxyons at different cen- 
trality; the experimental results are from [27], and the real 
lines are from our model. 



I m 
m 

° 80 



60 



40 



20 



1 Pb+Pb 


5.52TeV 




: LHC 







FIG. 3: The rapidity distribution of net baryons in central 

Fb+Pb collisions at LHC energies of y^sjvjv = 5.52TeV. The 
theoretical distribution is from our discussion. 



In Figure. 4, the star (*) is our prediction result of the 
mean rapidity loss < 6y > for Pb + Pb central collisions 
at LHC energies of = 5.52 TeV {yt = 8.68). In 

high-energy central nucleus-nucleus collisions, the baryon 
matter will be slowed down and will lose a few units of 
rapidity. The term "nuclear stopping power" was intro- 
duced in high-energy nucleus-nucleus collisions by Busza 
and Goldhabcr to refer to the degree of stopping an in- 
cident nucleon suffers when it impinges on the nuclear 
matter of another nucleus. Usually the mean rapidity 
loss < Sy >= 7/6— < y > was used to represent the nu- 
clear stopping power of nucleus -nucleus collisions. 

(4) It is shown in Fig. 5 that the dependence of per- 
centage ratios of the net-baryon production from the 
central gluon saturation region on incident energy at 
^/s^ = 17.3, 62.4 , 200 and 5520 GeV. It is found 
that the percentage ratio from the central gluon satu- 
ration region increases with the increase of the colliding 




FIG. 4: The dependence of mean rapidity loss < Sy > on 
beam rapidity yb- The A , •, ■, * are the calculated results 
by our discussion for ^/Inn = 17.3, 62.4, 200 and 5520 GeV, 
respectively. The real curve is the fit curve with < 5y >= 
1.548 • ln(j/6) + 0.036 
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FIG. 5: The dependence of the percentage of net-baryon from 
the central gluon saturation region on the colliding energies 
at y^sWN = 17.3, 62.4 , 200 and 5520 GeV. 



energies. The percentage ratio from the saturation re- 
gion rises rapidly from SPS y^sjviv = 17.3GeV to RHIC 
y/sWN = 62.4GeV, but slowly from RHIC ^/s^n = 200 
GeV to LHC ^/sJTn = 5520 GeV. It is found that the per- 
centage ratio is 43.17% at -^/sjvjv = 17.3GeV, but 56.31% 
at y/SNN = 5520 GeV. It seems that that the number of 
more than half of the net baryon at LHC comes from the 
central gluon saturation region. 

5The dependence of the effective quark mass rUq on 
incident projective rapidity in the CMS at -y/sjvAr = 17.3, 
62.4 , 200 and 5520 GeV is shown in Fig.6. It is found 
that the effective quark mass varies slowly with the 
varying of the incident energies, and rUq varies among 
0.24 to 0.26. 
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FIG. 6: The dependence of the effective quark mass on inci- 
dent projective rapidity in the CMS at y'sjviv ~ 17.3, 62.4 , 
200 and 5520 GeV. 

IV. SUMMARY AND CONCLUSION 

As discussed in Ref. [14,15], two distinct and symmet- 
ric peaks with respect to rapidity y occur at SPS energies 
and beyond in A -I- A colhsions. A geometrical scahng 
feature of peak position of net baryon rapidity distri- 
butions was propos ed to discuss about the net-baryon 
distribution [ij, [l^. The rapidity separation between 
the peaks increases with energy and decreases with the 
increasing mass number, A, reflecting larger baryon stop- 
ping for heavier nuclei, as was investigated phenomeno- 
logically in the Non-uniform Flow Model(NUFM) [lol - 
|l3|. In this work we show the geometrical scaling with 
gluon saturation rapidity limit, and also discuss the net- 
baryon rapidity distribution feature in the SPSRHIC and 
LHC. 

A saturation model for net-baryon distributions that 
successfully describes the net-baryon rapidity distribu- 
tions and their energy dependence is presented in this 
paper. The remarkable feature of geometric scaling pre- 
dicted by our discussion is reflected in the net-baryon ra- 
pidity distribution, providing a direct test of gluon sat- 



uration rapidity and x regions. The gluon saturation 
model is proposed by introducing a rapidity variable with 
gluon saturation region to define the gluon saturation re- 
gion of central rapidity region of centrally colliding heavy 
ions at ultra-relativistic energies. 

The gluon saturation features of central rapidity at 
SPS and RHIC can be investigated. It is found that 
the values of central rapidity of gluon saturation region 
increase with colliding energy. The detailed dependence 
of rapidity {y saturation) of central gluon saturation on 
colliding energy is also investigated in this paper. We also 
predict the net baryon rapidity distribution in central 
Pb-|-Pb collisions at LHC energies of ^Jsjsin = 
by gluon saturation model with geometric scaling. The 
gluon saturation region is larger than those of RHIC and 
SPS, and the separation of two symmetric peaks of net- 
baryon is much wider than that of SPS and RHIC. 

It is shown that gluon saturation feature is an impor- 
tant feature with the increasing of colliding energy. It 
seems that more than half of the produced net baryon 
numbers (56.31%) at LHC come from the central gluon 
saturation region, but the percentage ratio is 43.17% at 
SPS ./IWn = 17.3GeV. 

The dependencies of the percentage of net-baryon from 
the central gluon saturation region and stopping power 
on colliding energies are also studied in this paper. It 
is also found that the mean rapidity loss shows a linear 
dependence on \n{yb). From that we can predict the mean 
rapidity loss for future LHC experimental data. 



V. ACKNOWLEDGMENTS 

This work was supported by National Natural Science 
Foundation of China (10975091), Excellent Youth Foun- 
dation of Hubei Scientific Committee (2006ABB036)and 
Education Commission of Hubei Province of China 
(Z20081302).The authors is indebted to Prof. Lianshou 
Liu for his valuable discussions and very helpful sugges- 
tions. 



[1] Breitweg J et al. Eur. Phys. J. C, 1999, 7: 609 

[2] lancu E, McLerran M. Phys. Lett. B, 2001, 510: 145 [13] 

[3] Arsene I C et al. Phys. Lett. B, 2009, 677: 267 

[4] Bearden I G et al. Phys. Rev. Lett., 2004, 93: 102301 [14 

[5] Klay J L et al. Phys. Rev. Lett., 2002, 88: 102301 [15 

[6] Bearden I G et al. Phys. Rev. Lett., 2005, 94: 162301 [16 

[7] Barrette J et al. Phys. Rev. C, 2000, 62: 024901 [17 

[8] Ahle L et al. Phys. Rev. C, 1999, 60: 064901 [18 

[9] Appelshauser H et al. Phys. Rev. Lett., 1999, 82: 2471 

[10] FENG Shengqin, LIU Feng, LIU Lianshou. Phys. Rev. [19 

C, 2000, 63: 014901 [20^ 

[11] FENG Shengqin,YU AN Xianbao,SHI Yafei.Modern Phys [21 

LettA, 2006, 21: 663 [22 

[12] FENG Shengqin, XIONG Wei. Phys. Rev. C, 2008, 77: [23 



044906 

FENG Shengqin, YUAN Xianbao. Science in China Se- 
ries G, 2009, 52:198 

Mehtar-Tani Y et al. Phys. Rev. C, 2009, 80: 054905 

Mehtar-Tani Y et al. Phys. Rev. Lett., 2009, 102: 182301 

Gribov L V et al. Phys. Rep., 1983, 100: 1 

Mueller A H, Qiu J W. Nucl. Phys. B, 1986, 268: 427 

McLerran L, Venugopalan R. Phys. Rev. D, 1994, 50: 

2225 

Kharzeev D et al. Nucl. Phys. A, 2005, 747: 60 
Kharzeev D, Levin E. Phys. Lett. B, 2001, 523: 79 
Kharzeev D et al. Phys. Lett. B, 2004, 599: 23 
Baier R et al. 2006, Nucl. Phys. A., 764: 515 
Triantafyllopoulos D N. Nucl. Phys. B, 2003, 648: 293 



6 



[24] Kuracv E A ct al. Sov. Phys. JETP, 1977, 45: 199 [27] Debbe R et al. J. Phys.G. 2008, 35: 104004 

[25] Ya.Ya.Balitsky et al. Sov. J. Nucl. Phys, 1978, 28: 822 
[26] Lipatov L N. Sov. J. Nucl. Phys. 1976, 23: 338 



